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Abstract. 

We evaluate the Meteo-France global chemistry transport 3D model MOCAGE 
(MOdelc de Chimie Atmospherique a Grande Echelle) using the important 
set of aircraft measurements collected during the ICARRT/ITOP experiment. 
This experiment took place between US and Europe during summer 2004 
(July 15- August 15). Four aircraft were involved in this experiment provid- 
ing a wealth of chemical data in a large area including the North East of US 
and western Europe. The model outputs are compared to the following species 
of which concentration is measured by the aircraft: OH, H 2 O 2 , CO, NO, NO 2 , 
PAN, HNO 3 , isoprene, ethane, HCHO and O 3 . Moreover, to complete this 
evaluation at larger scale, we used also satellite data such as SCIAMACHY 
NO 2 and MOPITT CO. Interestingly, the comprehensive dataset allowed us 
to evaluate separately the model representation of emissions, transport and 
chemical processes. Using a daily emission source of biomass burning, we ob- 
tain a very good agreement for CO while the evaluation of NO 2 points out 
incertainties resulting from inaccurate ratio of emission factors of NO x /CO. 
Moreover, the chemical behavior of O 3 is satisfactory as discussed in the pa- 
per. 
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1. Introduction 

Modeling constitutes an essential complement to filling gaps in temporal and spatial 
coverage of in situ and remote sensing measurements. Numerous studies have shown the 
capacity of global chemistry transport models (CTM) to reproduce the main characteris- 
tics of the atmospheric chemical composition [Bey et al, 2001; Wang et al. , 1998; Horowitz 
et al. , 2003]. These successes have led to the use of CTM for chemical forecast [ Dufour 
et al. , 2005], satellite retrievals [ Palmer et al., 2001] and climate predictions [Dentener et 
al., 2006] 

The MOCAGE (MOdele de Chirnie Atmospherique a Grande Echcllc) model is a 3-D 
CTM which has been developped at the Centre National de Recherches Meterologiques 
(Meteo-France). Among current CTMs, MOCAGE presents several specificities. First, 
it is a multi-scale model, covering scales from regional with a high resolution down to 
0.1°x 0.1°to planetary with the resolution of 2°x 2°allowing the representation of small 
scale dynamic processes. Moreover, it combines tropospheric and stratospheric chemistry, 
which is important in e.g. ozone budget studies [ Riviere et al, 2005]. Finally, it uses 
a very detailed tropospheric chemical scheme RACM (Regional Atmospheric Chemistry 
Mechanim) [Stockwell et al, 1997]. This makes MOCAGE a very attractive model allow- 
ing it to cover a wide range of scientific objectives, from the chemical data assimilation 
[Pradier et al., 2006] to the modeling of tropospheric chemistry at a regional scale [ Dufour 
et al, 2003]. 
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During summer 2004, the European Intercontinental Transport of Ozone and Precursors 
(ITOP) was conducted as part of ICARTT [ Fehsenfeld , 2006]. Several scientific teams 
from Germany, France, the UK and the US carried out aircraft measurements of chemical 
species concentrations between N. America and Europe were involved in this experiment. 
The aim of this project was to better understand the mechanisms of pollution chemistry 
within up-lifted air exported from N. America to Europe. The campaign took place 
between July 15 and August 15, 2004. Strong wildfires over North America and Canada 
occurred with a high variability both in space and time and created rich pollution air 
plumes reaching the coasts of Europe and releasing numerous species in the atmosphere. 
The summer of 2004 was one of the strongest fire seasons on record for Alaska and Western 
Canada. The number of species simultaneously measured, the frequency of sampling, 
together with the large area considered, make this campaign a very useful data base for 
the evaluation of models. 

During this campaign, MOCAGE was used to provide chemical forecast and analysis 
to help defining the experimental aircraft flight plans. 

We present in section 2 the model set-up and the data used for this model evaluation 
(section 3). The results from the MOCAGE (Version 1.0) simulation in the period of the 
ITOP experiment are presented in section 4. In that section we compare modeled and 
observed concentrations of ozone and related trace gases using aircraft in-situ measure- 
ments and satellite data from the intruments MOPITT (Measurement Of Pollution In 
The Troposphere) [Drummond and Mand , 1996] and SCIAMACHY (SCanning Imaging 
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Absorption spectroMeter for Atmospheric CHartographY) [Bovensmann, H. et al, 1999]. 
Finally, we conclude in section 5. 

2. Model set-up 

The MOCAGE model is a global 3-D CTM providing numerical simulations of the in- 
teractions between dynamical, physical and chemical processes in the troposphere and 
lower stratosphere. MOCAGE uses a semi-lagrangian gridpoint model [ .Josse et al., 2004] 
to transport the species. For our simulation we used the global grid with a horizontal 
resolution of 2°x 2°. MOCAGE includes 47 hybrid (cr, P) levels from the surface up to 5 
hPa. The vertical resolution is 40 to 400 m in the boundary layer (7 levels) and about 800 
m in the lower stratosphere and in the vicinity of the tropopause. The chemical scheme 
used is RACMOBUS, which combines the stratospheric scheme REPROBUS [ Lefevre et 
al, 1994] and the tropospheric scheme RACM [Stockwell et al, 1997]. RACMOBUS in- 
cludes 119 individual species, among which 89 are prognostic variables, and considers 372 
chemical reactions. Convective processes are simulated with the scheme of Bechtold et al 
[2001], and turbulent diffusion is calculated with the scheme of Louis [1979]. MOCAGE 
also parameterizes emissions and dry deposition [Michou and Peuch, 2002; Michou et al, 
2005; Nho-Kim et al, 2004], and scavenging. MOCAGE uses the emission inventory from 
Dentener et al [2004, 2006] with a monthly or yearly resolution depending on the species. 
However, because of numerous wildfire events that occurred over N. America during the 
ICARTT/1TOP campaign, we modified the emissions of several species using the daily 
US emission inventory of Pfister et al. [2005]. The source emissions are just injected from 
the surface. For more details concerning the sources inventories used in the model see 
Table 1. The meteorological analyses of Meteo-France ARPEGE [ Courtier et al, 1991] 
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were used to initialize and constrain the dynamics of the model every 6 hours. The ver- 
tical velocity is calculated from the ARPEGE wind horizontal components by imposing 
the mass conservation law for each atmospheric column. 

Figure 1 shows the MOPITT CO concentrations binned at 0.5°x 0.5°at 500 hPa during 
the ITOP period (2004/07/15 to 2004/08/15), along with the corresponding horizontal 
wind held obtained from ARPEGE meteorological analyses. We see typical summertime 
features of the atmospheric circulation with a strong high above the western Atlantic and 
a low near the British Islands. Because of its long residence time, CO is known as a good 
tracer of polluted air masses and clearly shows here the pathway of exported pollution 
over the Atlantic, centered at about 50°N. Maximum CO concentrations are localized over 
Alaska and Canada due to the wildfire events during this period. The strong zonal how 
above Alaska and the following deep thalweg allow a rapid transport of pollutants to the 
N. E. US. The polluted plumes subsequently undergo a strong South West outflow. Our 
simulation started from a climatological initial held on June 1, 2004 at 12 UT. 

3. Rationale of the evaluation 
3.1. Chemical species 

We focused our evaluation of MOCAGE on the following species: OH, H 2 0 2 , CO, NO, 
N0 2 , PAN, HN0 3 , Isoprene, Ethane, formaldehyde (HCHO) and 0 3 . OH is the main 
oxidant for non-radical species in the atmosphere. Its concentration is of primary impor- 
tance for quantifying chemical processes in the troposphere, in particular the formation 
and destruction of 0 3 . The radical OH is formed by 0 3 photolysis as follow: 
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0 3 + hv 0 2 + 0( 1 D) 

0( 1 D) + H 2 0 -»• 20H 

Besides its fundamental role in tropospheric chemistry as principal OH source, the O 3 
molecule is known as being toxic to humans and vegetation because it oxidizes biological 
tissue. In the troposphere, O 3 can be produced by oxidation reactions between OH and 
some trace gas constituents. Moreover, H 2 O 2 is a reservoir for HO x (HO x = H+OH+HO 2 ) 
in the atmosphere. 

CO produced by incomplete combustion of hydrocarbons, plays a key role in tropo- 
spheric chemistry as ozone precursor. The chain mechanism [Jacob, 2000] leading to O 3 
production consists in the oxidation of CO by O 2 which is catalyzed by the HO x fam- 
ily and by NO x = NO + NO 2 . The lifetime of CO at mid-latitudes is about one month. 
Oxidation of hydrocarbons follows the same scheme as CO oxidation and is of primary im- 
portance in the regulation of tropospheric OH and O 3 concentrations. Ethane is released 
by industrial and combustion sources and is removed by OH oxidation with a lifetime of 
a few months. HCHO is a by-product of hydrocarbons oxidation. It can be scavenged 
by clouds but its lifetime is long enough to allow transport in remote region where its 
photolysis produces HO x radicals. The principal biogenic hydrocarbon contributing to 
O 3 formation is isoprene, an odorless compound which is a by-product of photosynthesis. 
Isoprene reacts extremely rapidly with OH, resulting in an atmospheric lifetime of less 
than one hour. NO x is produced by combustion processes, lightning and soil decomposi- 
tion. NO x concentrations control to a large part the chemical equilibriums involving O 3 


DRAFT 


May 31, 2006, 10:12pm 


DRAFT 



X - 8 


BOUSSEREZ ET AL.: MOCAGE EVALUATION 


production or destruction in the troposphere and thus play a key role in air quality. The 
principal sink of NO x is oxidation by HNO3. Because of its high solubility in water, HNO3 
is scavenged by precipitation in the troposphere. 

PAN is produced in the troposphere by photochemical oxidation of carbonyl compounds 
in the presence of NO x . In turn, PAN can regenerate NO x by thermal decomposition. 
The lifetime of PAN strongly depends on temperature, varying from 1 hour at 295 K to 
several months at 250 K. In the lower troposphere, NO x and PAN are near chemical equi- 
librium. However, in the middle and upper troposphere PAN is a reservoir for NO x . It 
can be transported over long distances and decomposed to release NO x far from its source. 

The ICARTT/ITOP campaign provided in situ measurements of all these species. In 
addition, satellite data allow global observations of pollution. 

3.2. The data 

The measurements were made aboard several aircraft: FAAM BAE 146 from the Insti- 
tute for Atmospheric Science (UK), DC-8 from the NASA Project Office (USA), P3 from 
the NOAA Aeronomy Laboratory (USA), and Falcon from the Institute of Atmospheric 
Physics (Germany) (see Table 2). MOZAIC (Measurement of OZone and water vapour 
by Airbus in-service airCraft) data [ Thouret et al, 1998; Nedelec et al, 2003] during the 
ITOP period are also used. We averaged the datasets over 1-min intervals and performed 
an on-line comparison to compare modeled and aircraft measurements. The observed 
and simulated concentrations obtained by this method were subsequently averaged over 
three regions (Fig. 2) : N. E. America, Europe and Atlantic ocean, referred to hereafter 
as domains 1, 2 and 3, respectively. The number of domains used depend also on the 
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availability of the aircraft data over the area. In order to evaluate our simulation at a 
larger scale, we also incorporated MOPITT CO retrievals and NO 2 tropospheric columns 
from SCIAMACHY. 

MOPITT is a nadir Infra-Red correlation radiometer onboard the NASA Terra Satel- 
lite. It has a horizontal resolution of 22 kmx22 km and provides a global coverage in 
about 3 days. We used the Level 2 V3 MOPITT datasets, which consist of retrieved 
CO mixing ratios for 7 vertical levels in the atmosphere (surface, 850 hPa, 700 hPa, 500 
hPa, 350 hPa, 250 hPa and 150 hPa). A detailed description of the MOPITT-CO re- 
trievals, based on radiance inversion, is given in [ Deeter et al, 2003]. The number of 
degrees of freedom in one profile, which characterizes the information content of the re- 
trieval, is typically less than 2 [ Deeter et al, 2004], For the comparisons with the model 
results, we applied the MOPITT averaging kernels to make comparable the different data. 

SCIAMACHY observes the upwelling radiation from the earth surface and the extrater- 
restrial solar radiance. It alternately measures in nadir and limb, covering the 220-2240 
nm spectral region range with a resolution of 0.25 nm in the UV, 0.4 nm in the visible 
and less than 0.4 nm in the NIR. The typical size of the nadir ground-pixel for NO 2 is 30 
km x 60 km. Its swath width is 960 km, providing a global coverage at the equator within 
6 days. Here we used the N0 2 tropospheric column product. The retrieval approach used 
for NO 2 nadir measurements is based on the Differential Optical Absorption Spectroscopy 
(DOAS) method. Details on the data analysis can be found in Richter et al., [2005]. 

Note that an on-line comparison was applied to compare MOPITT CO retrievals with 
MOCAGE CO simulations. In addition, comparisons of NO 2 columns between SCIA- 


DRAFT 


May 31, 2006, 10:12pm 


DRAFT 



X - 10 


BOUSSEREZ ET AL.: MOCAGE EVALUATION 


MACHY and MOCAGE were carried out using an off-line interpolation. The MOCAGE 
tropospheric columns of N0 2 are calculated from the surface to the tropopause as deter- 
mined with a PV = 2. 

4. Results and discussions 

4.1. Hydroxyl Radical and Hydrogen Peroxide 

Figures 3 and 4 present the vertical profiles of OH and H 2 0 2 respectively. Over do- 
main 1 the model slightly overestimates OH concentrations by 0. 1-0.2 pptv in the lower- 
troposphere. OH concentrations are very sensitive to anthropogenic emissions close to 
the surface. The model divergence is thus probably linked to inaccuracies in our emis- 
sion inventory. The good agreement with observations at mid-troposphere attests that 
the tropospheric chemical equilibrium is well simulated by MOCAGE. Above 8 km, the 
underestimation of 0. 1-0.2 pptv may be due to an underestimation of cloud albedo pho- 
tochemical effects [ Brune et al., 1998] or too low H 2 0 concentrations. Simulated H 2 0 2 
concentrations are slightly lower (bias of 500 pptv) than observed concentrations in the low 
and the middle troposphere of domain 1 but remain within the range of observations. The 
model reproduces well the decrease of H 2 0 2 with increasing altitude due to the airmass 
drying out. Lower H 2 0 2 concentrations from the model between 0 and 5 km are consistent 
with higher OH concentrations that enhance the reaction: H 2 0 2 + OH — > H0 2 + H 2 0. 
Over domain 3, the model overestimates the observed OH concentrations by about 0.36 
pptv throughout the troposphere. The bias in lower-troposphere may result from an un- 
derestimation of the marine boundary layer cloud cover in ARPEGE [ Mathieu et al, 2003] 
leading to a higher rate of ozone photolysis. More globally, the OH tropospheric bias may 
be linked to photochemical effects of aerosols not included in MOCAGE [ Duncan et al., 
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1997]: (1) aerosol scattering and absorption of UV radiation, and (2) reactive uptake of 
H0 2 , N0 2 and N0 3 . For example, Martin et al. [2003] showed that aerosols may decrease 
boundary layer OH concentrations by a factor 2 over northern Europe during summer 
and up to a factor 4 over biomass burning regions of North India. During the ITOP 
period, long range transport (LRT) events were precisely associated with high aerosol 
concentrations in pollution plumes [ Real et al., 2006]. Results for H 2 0 2 over domain 3 
show an underestimation of about 1000 pptv consistent with the higher oxidizing power 
of the atmosphere in the model. 

4.2. Carbon Monoxide 

Figure 5 shows the comparisons between MOCAGE and MOPITT CO at 500 hPa for 
the period July 15 - August 15, 2004. The model reproduces very well the observed distri- 
bution. The main pollution pathway spreads from N.E. America to N.W. Europe, centered 
around 50°N. The correlation coefficient (r 2 ) between the model and the observations at 
this level is approximately 0.98. The histograms (Fig. 5 c) are very comparable though 
the model tends to slightly overestimate the high and the low values. We calculated the 
figure of merit defined as the ratio of the surface of the histograms intersection and the 
surface of the histograms union. Closer to 1 closer the distributions are. The calculated 
value is 0.66. Figure 6 shows MOPITT and MOCAGE CO at 500 hPa for the period 
between July 24 and July 26, 2004. This case shows that MOCAGE allows an accurate 
representation of relative small-scale structures seen by MOPITT, in particular over the 
Atlantic ocean and N.W. Europe. 
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Figure 7 displays the aircraft and simulated vertical profiles of CO. Observed profiles 
are similar over domain 1 and 2 and fairly well reproduced by the model. Maximum 
values are found in the lower-troposphere due to anthropogenic emissions. Convection 
processes over continental areas during the summer result in a well mixed middle and 
upper troposphere. There is a weak CO concentration overestimation of about 10-20 
ppbv above 4 km globally remaining within the range of observations. Results for domain 
3 show approximately the same bias throughout the troposphere. This bias which is a 
weak overestimation of modeled CO background is probably due to an overestimation of 
CO emissions in the model. CO comparisons over domain 3 establish the capability of 
the model to reproduce accurately the transport processes. Indeed, the two peaks at 6 
km and 9 km as well as the high concentration variability linked to LRT events are well 
reproduced by the model. 

4.3. Hydrocarbons and Isoprene 

Comparisons of vertical profiles of ethane are shown in Figure 8. For domain 1, there is 
a good agreement between simulated and observed concentrations above 2 km. The lower 
troposphere concentrations are underestimated by about 300 pptv but remain within the 
range of observations. Comparisons for domain 2 show the same behavior as domain 1. 
ffowever, very few data are available below 8 km altitude. 

The global emission inventory for ethane varies in the litterature between 8 and 36 
Tg(C).y _1 . The inventory used in MOCAGE calculates an annual emission rate of 8 
Tg(C).y _1 . ffowever, Xiao [2005] found an estimate of 13.5 Tg(C).y _1 . Thus, the bias 
between measured and modeled concentrations in the continental Planetary Boundary 
Layer (PBL) might be due to an underestimation of ethane emissions in the model. In 


DRAFT 


May 31, 2006, 10:12pm 


DRAFT 



BOUSSEREZ ET AL.: MOCAGE EVALUATION 


X - 13 


addition, the higher variability of the observations relative to the simulated values near 
the ground shows the strong dependence of ethane to emission resolution. For domain 3, 
the simulated ethane profile is consistent with the observed profile except between 9 and 
10 km. The very few data available at those altitudes possibly cause a representativeness 
error. However, the model reproduces well the two peaks at 3 and 6 km as well as the 
high variability of the concentrations. 

The HCHO comparisons between model results and aircraft observations are presented 
in Figure 9. Over domain 1, the simulated profile is consistent with the observed profile 
despite a slight overestimation of about 100 pptv in the middle and upper troposphere. 
Over domain 3, the shape of the HCHO profile is quite well reproduced by the model. 
There is a non-systematic weak underestimation outside the PBL probably resulting from 
too high OH concentrations that shorten the HCHO lifetime by oxidation reactions. Note 
that oxidation by OH radicals leads to a CO production which can be linked to the CO 
bias discussed previously. The high variability of observed values compared with the low 
variability of simulated values reflects these lifetime differences. Indeed, a longer lifetime 
allows a greater influence of LRT events on concentration profiles. Besides, we note for 
the observations a nearly constant value with altitude of the order 300 pptv above 5 km, 
which has been identified in previous studies as the signature of air masses originating from 
North America [ Helan et ai, 2003] Figure 10 presents the results of the comparison for 
isoprene over domain 1. Model concentrations are in quite good agreement with observed 
concentrations except near the surface where the strong dependence of isoprene emissions 
to soil nature and vegetation induces a high variability. 
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4.4. Nitrogen species 

Figure 11 compares the NO 2 tropospheric column retrieved by SCIAMACHY and that 
simulated by MOCAGE for the period between July 15 and August 15, 2004. Globally, the 
model overestimates NO 2 column concentrations by about 3.10 14 molec.cm -2 . Figure 11c 
shows that the most pronounced biases are localized over low NO 2 concentration regions 
i.e. over oceans. However, the shape of histograms from the model and the observations 
are quite similar with a figure of merit of about 0.67. The NO 2 global held is quite well 
reproduced by the model except for N.W. America under strong biomass burning. 

Turquety et al., [2006] recently estimated that 17 % of the total area burned in Alaska 
and Canada during the summer 2004 was located in peatlands. In addition, Bertschi et al. 
[2003] studied the emissions from residual smoldering combustion (RSC) and concluded 
that correcting emission estimates by including RSC could decrease the NO emission 
factor (EF) by a factor 2. We found that our EF(NO x )/EF(CO) ratio is more than 
twice that recommended by Bertschi et al. [2003], which could explain the discrepancies 
observed between SCIAMACHY measurements and MOCAGE simulation over wildfire 
areas. Finally, the NO 2 concentrations simulated over the Atlantic are very high com- 
pared to SCIAMACHY measurements. Note that Richter and Burrows [2002] have shown 
that most of the uncertainties in the derived tropospheric columns using DOAS method 
result from a contamination by clouds and light path, leading in some cases to an un- 
derestimation by up to 40 %. The difference in calculation of tropospheric columns may 
also explain the bias between MOCAGE and SCIAMACHY NO 2 . Moreover, the North 
Atlantic region is known as being very cloudy during summer [ Mathieu et al., 2003] and 
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this could considerably degrade the accuracy/precision of the retrievals. 

A comparison of simulated and observed NO profiles is shown in Figure 12. Over do- 
main 1 , the model is in relatively good agreement with the observations for the altitude 
range 0 to 5 km, while it underestimates the NO concentrations in the upper troposphere 
by up to 400 pptv. The measured NO profile is typical of lightning activity [ Decaria et 
al . , 2005; Barthe et al . , 2006]. And indeed, strong lightning activity occurred over domain 
1 between July 15 and August 15, 2004 (Fig. 14). Thus the upper troposphere NO bias is 
likely due to the fact that MOCAGE does not parameterize NO x production by lightning. 
Over domain 2 and 3, the simulated NO concentrations are consistent with measurements 
in the low and mid-troposphere. The slight divergence observed in the upper troposphere 
is probably also linked to NO production by lightning but is smaller than over domain 1 . 

Figure 13 shows the simulated and observed NO 2 concentration profiles. Over domain 1, 
NO 2 concentrations are fairly well reproduced by the model. The weak underestimation in 
the upper troposphere is less pronounced than for NO because of the very low air density 
and cold temperatures at this altitude that considerably reduce the speed of the kinetics 
of the reaction: NO + O 3 — > NO 2 + O 2 . Over domain 3, the model underestimates the NO 2 
concentrations by about 200 pptv throughout the troposphere. These discrepancies might 
be explained by the higher OH concentrations in the model resulting in a shorter NO 2 
lifetime. Besides, the high variability of the measurements at mid-troposphere compared 
to the simulated concentrations also seems to indicate a significant influence of LRT events 
on the measured profiles. Over domain 2, there is a relatively good agreement between 
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observed and simulated NO 2 concentrations in the lower troposphere except close to the 
surface where the high variability of measurements makes the comparison difficult. The 
systematic underestimation of about 200 pptv above 3 km may be due to intercontinental 
transport, based on the similarity with the bias observed over domain 3. 

Concerning PAN concentrations (Fig. 15), the model overestimates observed concentra- 
tions by about 500 pptv throughout the troposphere over domain 1. These discrepancies 
remain within the range of observations in the PBL, where the concentrations are essen- 
tially influenced by surface NO x emissions. Above 2 km, the bias probably results from 
too high NO x emission rates over the wildfire areas as pointed out previously, and the 
subsequent excess of PAN production. Over domain 3, the overestimation is quite similar 
to that observed over domain 1, except in the lower troposphere where PAN decomposi- 
tion takes place. The model reproduces both the peak at 6 km and the high variability 
in mid-troposphere linked to LRT events. Over domain 2, simulated PAN concentrations 
are consistent with observations at 1 km. However, above 1 km, a systematic bias quite 
similar to that found over domains 1 and 3 is observed, suggesting the influence of inter- 
continental transport of pollution. Note that the discrepancies observed near the ground 
can be due to a representativeness error owing to the very few data available. 

For HNO3 (Fig. 16), simulated concentrations are similar to observed concentrations 
over domain 1. As expected, the HNO3 profile is correlated with the NO2 profile as HNO3 
is produced in the troposphere by oxidation of NO 2 . For domain 3, observed concen- 
trations are overestimated by the model in the lower troposphere (from 0 to 3 km) up 
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to a factor 8 close to the surface. These strong discrepancies have already been men- 
tioned in previous model evaluations [ Horowitz et al . , 2003] and might be linked to the 
wet deposition parameterization and inaccurate rain fluxes in the meteorological analyses. 
In addition, the underestimation of the marine boundary layer cloud cover in ARPEGE 
analyses may also induce a slower wet deposition near the surface. Moreover, the high 
OH level simulated in the model and the presence of relatively high NO2 concentrations 
could lead to significant production of HNO3. Above 2 km, observed HNO3 concentrations 
are relatively well reproduced by the model except in the upper troposphere where cloud 
scavenging is probably underestimated. 

In order to better understand the chemical processes involved during the transport 
of polluted air masses, it is useful to consider the NO y family (NO y = NO + NO2 + 
PAN + HNO3). Figure 17 shows NO y profiles comparisons between observations and 
simulations. Over domain 1, the shape of the profile is relatively well reproduced. The 
model overestimates observed concentrations by about 500 pptv in the low and mid- 
troposphere while there is a slight underestimation in the upper troposphere probably 
due to NO x production by lightning. Over domain 3, there is approximately the same 
NO y overestimation as for domain 1 and we observe a fairly good correlation between 
the observed and modeled profiles except in the PBL and upper troposphere. Because 
nitrogen oxides from combustion processes and lightning are mainly emitted as NO x and 
converted into NO y species afterwards, the [NO x ]/[NO y ] ratio can be used as a measure 
of the age of emissions [ Helan et al, 2003]. 
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Figure 18 presents the NO y partitioning over domain 1 and 3 for both the simulation 
and the observations. Over domain 1, the global partitioning is quite similar in the model 
results and the observations. In particular, in the lower troposphere, the contribution of 
PAN rapidly increases with the altitude while HNO 3 contribution decreases. We also ob- 
serve that the modeled [NO x ]/[NO y ] ratio decreases with increasing altitude, whereas the 
observed [NO x ]/[NO y ] ratio increases in the upper troposphere where the lightning NO x 
production takes place. Over domain 3, the simulated NO y distribution is approximately 
the same as domain 1. Here the weak [PAN]/[NO x ] ratio near the surface is likely clue 
to the thermal decomposition of PAN. The contribution of NO 2 over domain 3 for the 
observations is higher than model results, which probably reflects differences in air masses 
reactivity resulting from differences in OH concentrations. 

4.5. Ozone 

Figure 19 displays O 3 profiles from MOCAGE and from aircraft measurements. The 
model overestimates O 3 concentrations by about 30 ppbv in the lower troposphere (below 
2 km) of domain 1. Several models found a comparable overestimation over the N.E. US, 
partly resulting from recent changes in isoprene and anthropogenic emissions [ Fiore et al, 
2005] and from the effect of the horizontal resolution. The bias decreases with increasing 
altitude and observed and simulated concentrations are in good agreement above 5 km. 
Over domain 3, the shape of the observed profile is fairly well reproduced despite a mean 
overestimation of about 18 ppbv throughout the troposphere. Chandra et al. [2004] found 
that over Northern oceans the contribution to the Tropospheric Column Ozone due to 
stratosphere-troposphere exchange is about 50-55 %. Thus, the relatively good repre- 
sentation of O 3 concentrations in the mid- and upper troposphere might be the sign of 
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efficient coupling of tropospheric and stratospheric chemistry in MOCAGE. Over domain 
2, the shape of the observed profile is quite well reproduced by the model despite a slight 
overestimation by about 20 ppbv in the lower troposphere. 

Several studies have shown relatively robust correlations between O3 and NO y [Helan 
et al, 2003]. We thus calculated A0 3 / ANO y ratios, where A represents the difference 
between modeled and observed concentrations. Over domain 1, we found a AO3 / ANO y 
ratio of about 10 below 2 km and about 33 between 2 and 3 km altitude. These values are 
in good agreement with those reported in the litterature for this area [ Helan et al., 2003]. 
Over domain 3, we established a AO3 / ANO y ratio of about 50 in the mid-troposphere, 
which corresponds to the values reported in Helan et al. [2001] for marine environment. 
Simulated 0 3 concentrations over domain 2 are relatively consistent with the observa- 
tions despite a slight overestimation of about 10 ppbv in the mid- and low troposphere. 
The significant underestimation occurring at 12 km is probably linked to a too elevated 
tropopause height in the model. 

Therefore, onr results show that NO y overestimations explain quite well the bias ob- 
served in O3 concentrations and suggest that the model reproduces relatively well the 
chemical equilibriums between O3 and the NO y family. Note that among the chemical 
processes leading to an ozone enhancement in pollution plumes during LRT events, the 
PAN decomposition is identified as a possibly dominant source as pointed out in Hudman 
et al. [2004], Thus, the slight O3 overestimation over the Atlantic domain may be linked 
to the excess of PAN in the model. 
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5. Conclusion 

This paper provides a comprehensive evaluation of the global 3D MOCAGE CTM us- 
ing aircraft measurements and satellite observations available during the ICARTT /ITOP 
experiment. Despite some significant discrepancies, the model is capable of representing 
the general features of ozone and related trace gas profiles over distinct regions. To per- 
form this evaluation, we defined 3 domains for comparing the model results to the aircraft 
observations: over N.E. US, Europe and the Atlantic corresponding to a pollution remote 
region, the area affected by this pollution and a region between where the pollution is 
transported, respectively. In addition, simulated CO and NO 2 are compared to satellite 
instruments MOPITT and SCIAMACHY data, respectively. We obtain the following re- 
sults: 

1. OH concentrations are relatively well simulated above the boundary layer of N.E. US. 
The overestimation of about 0.1 pptv observed at lower troposphere over this domain could 
reflect inaccuracies in our emission inventory. The model systematically overestimates 
observed OH concentrations by about 0.2-0. 4 pptv over the Atlantic domain. This bias 
may result from both an underestimation of the Atlantic cloud cover by the meteorological 
ARPEGE analyses and the lack of photochemical effects of aerosols in MOCAGE. 

2. The vertical behavior of H 2 O 2 over N.E. US is correctly represented by the model. 
The H 2 O 2 overestimation observed over the Atlantic domain is consistent with the OH 
bias. 

3. The model reproduces very well the CO fields. However, there is a slight system- 
atic overestimation of about 10-20 ppbv in the free troposphere which could reflect a 


DRAFT 


May 31, 2006, 10:12pm 


DRAFT 



BOUSSEREZ ET AL.: MOCAGE EVALUATION 


X - 21 


global overestimation of CO emissions. Comparisons with MOPITT CO suggest that the 
transport processes are particularly well simulated by MOCAGE. 

4. The model reproduces quite well the ethane profiles but an underestimation of about 
500 pptv is observed within the continental PBL that probably results from too weak 
emissions in our inventory. 

5. Comparisons for HCHO and isoprene show a relatively good agreement between ob- 
servations and simulation. Comparisons with aircraft measurements and SCIAMACHY 
NO 2 suggest that we overestimate NO x anthropogenic emissions over N.E. US and Europe. 
Above the boundary layer over N.E. US, NO x concentrations are well simulated except 
in the upper troposphere where NO x production by lightning is not yet parameterized 
in this version of MOCAGE. The underestimation of 100-200 pptv over the Atlantic and 
Europe domains is probably due to the excess of OH in the model. However, over the 
Atlantic the discrepancies remain within the range of observations. In the opposite, NO 2 
tropospheric columns from SCIAMACHY are overestimated by the model, which could 
reflect a large impact of the Atlantic cloud cover on the NO 2 retrieval errors. In addi- 
tion, comparisons with NO 2 from SCIAMACHY pointed out that the daily US emission 
inventory overestimates the NO x emissions over wildfire areas. This explains the near 
systematic bias in PAN concentrations of about 500 pptv in the model. HNO 3 profiles are 
relatively well reproduced except in the Atlantic marine boundary layer where the model 
probably underestimates wet scavenging processes. 

6 . O 3 profiles show an overestimation of about 30 ppbv in the lower troposphere of the 
N.E. US domain, which could result from recent changes in trace gas precursors emissions 
over the US. Over the Atlantic, the overestimation is about 10-30 ppbv. The simulated 
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ozone profile over Europe fits well with the observed profile. We note that globally, the 
ozone bias are well correlated with the NO y bias. 

Therefore, this study suggests that improvements are needed in our global and daily 
emission inventory, especially for NO x species, to better simulate the ozone and related 
trace gas concentrations. Further work will be done soon using a new version of MOCAGE 
with a lightning NO x parameterization and new emission factors. 
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Figure 1. MOPITT CO (ppbv) binned at 2 °x 2°and ARPEGE horizontal wind at 500 hPa 
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averaged from July 15 to August 15, 2004. 



Figure 2. Regions used to aggregate aircraft observations for the purpose of model evaluation. 
The different boxes represent 1: N.E. US region; 2: Europe region; 3: the Atlantic ocean region. 
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Figure 3. Comparisons of observed and simulated vertical profiles of OH concentrations (pptv). 
The open squares are mean observed values (with horizontal bars for standard deviation), and 
the open triangles and solid lines are median observed values. Open diamond are mean simulated 
values and cross and dotted lines are median simulated values (with horizontal bars for standard 
deviations). Values on the right are the number of aircraft data averaged for each level, a): 
domain 1; b): domain 3. 




Figure 4. Same as figure 3 (pptv) but for H 2 O 2 . 
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CO(ppbv) 

Figure 5. Comparison between MOPITT and MOCAGE CO (ppbv) during ITOP at 500 


hPa. a) MOPITT CO binned at 0.5 °x 0.5 °; b) MOCAGE CO; c) corresponding histograms of 


MOCAGE (dashed line) and MOPITT CO (solid line). 
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Figure 6. Comparison between MOPITT and MOCAGE CO (ppbv) at 500 hPa for the period 
between July 24 and July 26, 2004. a) MOPITT CO binned at 0.5 °x 0.5 °; b) MOCAGE CO. 
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Figure 7. Same as figure 3 but for CO (ppbv) and for a) domain 1; b) domain 3; c) domain 2. 
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Figure 8. Same as figure 7 but for ethane (pptv). 
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Figure 9. Same as figure 3 but for formaldehyde (HCHO) (pptv). 
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Figure 10. Same as figure 3 but for isoprene (ppbv) and only over domain 1. 
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Figure 11. Comparison between SCIAMACHY and MOCAGE N0 2 for the period between 
July 15 and August 15, 2004. The color code represents the logarithm of N0 2 concentration 
fjn^il^p.pmf 2 ). Top: N0 2 troposph<j|'^c 9 |imrj](§o%oniL(^GL^^(IACHY; Middle: N0 2 trgip^s^hprip 
columns from MOCAGE. Bottom: corresponding histogram of MOCAGE N0 2 (dashed line) and 
SCIAMACHY N0 2 (solid line). 
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Figure 14. Lightning activity over the US from Jnly 15 to August 15, 2004. The black box 
delimits the domain 1. Colored pixels indicate the number of cloud-to-ground lightnings averaged 
over a 0.5 °x 0.5 “grid, (from the National Lightning Detection Network) 
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Figure 15. Same as figure 7 but for PAN (pptv). 
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Figure 16. Same as figure 3 but for HN0 3 (pptv). 
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Figure 17. Same as figure 3 but for NO y (pptv). 
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Figure 18. Repartition of the different species for NO y (%) for the MOCAGE model and the 
aircraft measurements: a) model over domain 1; b) model over domain 3; c) observations over 
domain 1; d) observations over domain 3. 
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Table 1. MOCAGE emissions inventories used for the ITOP simulation. 


“Species 

Biomass burning 

Other sources 

CO 

(b) 

(c) 

NO x 

(b) 

(c) 

TOL 

(c) 

(b) 

ch 4 

(b) 

(c) 

hc 5 

(b) 

(b) 

ALD 

(b) 

(b) 

HCHO 

(b) 

(b) 

hc 3 

(b) 

(b) 

ETH 

(b) 

(b) 

ETE 

(b) 

(b) 

KET 

(b) 

(b) 

OLI 

(b) 

(b) 

OLT 

(b) 

(b) 


(a) see Stockwell et al., [1997] for the acronyms 

(b) Daily emissions over the US from Pfister et al., [2005] 

(c) from Dentener et al. [2004, 2006] 


Table 2. Chemical species measured by the different aircraft 


BAE-146 

P3 

DC-8 

Falcon Airbus (MOZAIC) 

CO, 0 3 

CO, 0 3 

CO, 0 3 

CO, 0 3 CO, 0 3 

NO, N0 2 

NO, N0 2 

NO, N0 2 

NO 

PAN 

PAN, HN0 3 

PAN, HNO s , H 2 0 2 , OH 


ETH, ISO, HCHO 

ETH, ISO, HCHO 

ETH, ISO, HCHO 
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